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Abstract. We show that an initially unpolarized electron flow acquires spin 
polarization after passing through a lateral barrier in two-dimensional (2D) system 
with spin-orbit interaction (SOI) even if the current is directed normally to the barrier. 
The generated spin current depends on the distance from the barrier. It oscillates with 
the distance in the vicinity of the barrier and asymptotically reaches a constant value. 
The most efficient generation of the spin current (with polarization above 50%) occurs, 
when the Fermi energy is near the potential barrier maximum. Since the spin current 
in SOI medium is not unambiguously defined we propose to pass this current from the 
SOI region into a contacting region without SOI and show, that the spin polarization 
loss under such transmission can be negligible. 
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1. Introduction 

Generation and manipulation of spin-polarized carriers in semiconductor structures 
solely by electric methods is a key problem of spintronics [H [21 [3] . One of widely studied 
approaches to attain this goal is based on using spin-orbit interaction (SOI). The SOI is 
known to produce the spin polarization of electron current in layered tunnel structures. 
The effect is caused by the Rashaba SOI at the barrier boundaries of asymmetric 
structures or by the Dresselhause SOI in the barrier bulk |5l El 13 El E] . A general 
property of such structures is the absence of the spin polarization in the case where the 
current is directed normally to the barrier. In other words, a current component along 
a barrier should be created to get a spin current. This limits the capability of these 
structures to generate spin currents. Recently it has been found that two-dimensional 
(2D) structures with a lateral barrier are free of this restriction |10j. The electron 
current passing through the barrier acquires spin polarization, which exceeds 50% even 
if the current is directed normally to the barrier. However, in the studied case the SOI 
exists only inside the barrier. Such structures seem to be hardly realizable, since it is 
problematic to localize the Rashba SOI within the lateral barrier, especially if the latter 
is created by gate electrodes. 

In the present work the research of Ref. |10j is generalized to the case when SOI 
exists everywhere: in the potential barrier and the surrounding electron gas. We find 
that the high spin polarization can also be achieved in such structures. However, in this 
case two important questions arise concerning the definition of the spin current and the 
existence of equilibrium spin currents in 2D electron gas with SOI [11]. They provoked 
recently a wide discussion [121 [131 [HI [IS]. As regards the existence of equilibrium 
spin currents, this question is not essential for the barrier structures considered here 
for the following reason. The equilibrium spin current is known to be generated only 
within a narrow energy layer —Ego < E < 0, where E is electron energy. Ego is a 
characteristic energy of SOI [11]. If the barrier height U considerably exceeds Ego the 
barrier transparency for electrons in this energy layer is negligibly small. 

The problem of the spin current definition can be overcome by calculating an 
observable physical quantity, which is well defined and closely related to the spin flow in 
SOI medium. This could be a spin current in a normal 2D electron gas (without SOI), 
which is brought to the contact with the SOI structure under consideration. In other 
words, it is reasonable to explore a structure in which the spin current generated in the 
SOI region passes into a normal region where spin current is unambiguously defined. 
One can say that this region is designed to simulate, at least partially, a measuring 
process. With this in mind we study the spin current transformation when electrons 
pass through a contact between SOI and normal 2D regions, and find conditions under 
which this transformation occurs practically without loss of spin polarization. 

Finally we have found, that the barrier in the 2D electron system with the SOI 
allows one to generate the electron current with spin polarization exceeding 50% and 
this spin current can be transferred into a normal 2D electron gas with minimal loss. 
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Figure 1. Complex band structure of 2D electron gas with SOI. Energy branches 
1±)2± are shown as functions of k'^ and fc". Branch 3 is defined along real energy 
trajectories in a complex plane (fc^,fc"); only one branch located in the quadrant 
(fc^, /c" > 0) is represented. Panels a and b correspond to cases ky < a and ky > a. 



2. Basic wave functions and energy spectrum 

The structure to be studied here is a sheet of 2D electron gas with Rashba SOI separated 
by a potential barrier into two semiplanes (reservoirs) between which a small voltage 
V is applied. We are going to find electron and spin currents through the barrier, but 
begin with a discussion of wave functions for the whole system. The system is described 
by the Hamiltonian: 

H = ^^-^^ + '^{Pya^-p^ay) + U{x), (1) 

where components of electron momentum, a is the SOI parameter, a^^y the 

Pauli matrices and U{x) the barrier potential. We consider here a rectangular barrier 
of high U and width d\ U{x) = U at < x < d and U{x) = at a; < and d < x. The 
effective mass is supposed to be independent of the coordinates. 

Wave functions in the barrier and reservoirs are presented in the form of a linear 
combination of basic eigenfunctions of homogeneous 2D electron gas with SOI 

= E 

s' 

where s is the spin index, k = {kxs,ky) is the wave vector, k = {—kxs,ky). The 
wave vector component k^s is different for the barrier and the reservoirs. In addition 
it depends on the spin. In contrast the component ky is a conserved quantity and 
hence it is the same for all regions in a given state. The eigenfunctions and energy 
spectrum of homogeneous electron gas were studied in detail in Ref. [TD]. The main 
results which will be used below, are the following. Since the considered system is 
not translationally invariant in the x direction, the wave vector component k^s can be 
complex: k^ = k'^ + ik". In contrast, ky is always real. The total spectrum includes 
three spin-split branches (see Fig. 1). 



xAK 
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1) The first branch 1± corresponds to propagating states {k'^ = 0) with energy: 

Ck,s = -a'+(a + sy^27^ , (3) 
and spin function: 

where Ck,s = '^E]^,s'^/^'^, -S'k.s is the electron energy, a = ma/h? is the characteristic wave 
vector of the SOI, Cso = o,^ corresponds to the characteristic energy. Ego = h^a^ /2m. 

2) The second branch 2± exists when ky ^ 0, in the energy gap between branches 
1+ and 1_. These states decay monotonously with x and hence k'^ = 0. The energy and 
spin functions are defined by equations ([3]) and (jl]), where k'^ must be replaced by ik". 

3) The third branch lies below two above considered branches, Ck,s < -a^. It IS 
defined for real energy trajectories in the complex plain {k'^, k'!^): 

k'X + a\kl + k'^'-k':')-a' = 0. (5) 

The energy and spin functions for this branch are 

Ck,. = -a'-^, (6) 



a 



h — h" -V- ih' 

Note that at any given energy and ky, there are 4 eigenstates. In the case of first and 
second branches, the different eigenstates correspond to different signs of s and k'^ or 
k'l- For the third branch the eigenstates differ by signs (±) of k'^ and k'^. 



3. Spin polarized current through a barrier 

We now turn to the calculation of electron and spin currents flowing normally to the 
barrier. For simplicity suppose that the 2D electron reservoirs to the left and right 
of the barrier are equipotential and the potential difference V is small as compared 
to all characteristic energies of the system. The electron states contributing to the 
current are located in the energy interval of eV width near the Fermi energy Ep. In 
the {kx, ky) space, they occupy two semirings corresponding to electrons with opposite 
spins (flg. 2). The currents are determined by the summation of partial currents over 
these states [Ml [IT] . Using variables C and ky one flnds: 

w> X eV ^ kF,s f^^'^ dky 

J^^f) = ^2^ / ', , / / ^ j{CF,ky,s), (8) 

8^' V VCf + J~k,^s ^k%^ - kl 

where Cf = ^mEp/h^, kp^s is deflned by the equation: 

Cf = -a^ + (a + skp^sf, (9) 
and jXCf, ky, s) is the partial current in the eigenstate \(p, ky, s). 
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Figure 2. Energy diagram of the barrier structure. Lines 1±,2±,3 represent the 
spectrum branches described in the text. In the inset: full semirings in the {kx,ky) 
space which are occupied by electrons contributing to the current. 



The current J^Cf, ky, s) is calculated using the wave functions defined in equation ([2]) 
as a linear combination of basic eigenf unctions, the set of four eigenfunctions being 
different for the barrier and reservoirs as well as the spectrum there. The selection of 
basic eigenfunctions from all three space regions and all spectrum branches to form the 
total wave function corresponding to a given energy ( and transverse momentum ky is 
an intricate problem. Its solution is summarized in the diagram shown in Fig. 3. There 
are 12 regions on the plane {(, ky). The regions are bounded by four curves, 1-4, which 
are determined by the equations: 

ky= y/C + a^Ta, (10) 
ky = ^/C-u + a'^^a, (11) 
where u = 2mU/fi?. 

In each region a specified set of four eigenfunctions to be used in forming the total 
wave function is pointed. The list of these regions and corresponding eigenfunction sets 
for the reservoirs and the barrier are the following: 

1) For the reservoirs: Ri is the region without propagating states; R2 contains 
two waves of spectrum branch 1_, which are incident on the barrier, and two waves of 
branch 1_, which are refiected. For brevity we depict this schematically as follows: R2 
- (^ 1_, ^ 1_). The arrows designate the right- and left-moving waves, the number of 
arrows specifies the number of waves, and the figures behind them indicate the spectrum 
branches they belong to in accordance with Fig. 1. Using these notation, other regions 
are imaged as: R^ - {-^ 1-,^ 1-,^ 2_); Ra, - (^ 1_,^ l+i^ 1-;^ 1+); R5 ^ 
(-> 1_,^1_,^2+). 

2) For the barrier: Bi contains 4 modes of branch 3; i?2 ^ 4 modes of branch 1_; 
i?3 - 2 modes of 1_ and 2 modes of 2_; ,84-4 modes of 2_; -B5 - 4 modes of 1+; - 
2 modes of 1+ and 2 modes of 2+; B-j - A modes of 2+. 
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Figure 3. The distribution of the basic eigenfunctions in the plane of parameters 
{C,ky) for the reservoirs (R) and the barrier (B). The inner borders (hnes 1,2,3) 
determine the regions with different sets of four eigenfunctions. Thick hnes 4 and 
5 are external borders outside of which no propagating states exist. One more border, 
closing the region of the states accessible for electrons, is the Fermi energy (p. Dashed 
line at Ci corresponds to an example considered in the text. 



Fig. 3 helps to find the eigenfunction sets forming the total wave function for a given 
energy (i. It is needed to draw a vertical line C = Ci- The regions, which it crosses, 
show the eigenfunction sets according to the above list. If this line crosses more than one 
region, the integration interval in equation ([8]) is to be divided into parts corresponding 
to its intersection points with internal lines. 

As an example, let us describe the tunneling of electrons with energy E < U — Ego- 
If an electron falls on the barrier from the left reservoir in the state 5^, the wave 

function in this reservoir, x < 0, is 

\'^k2,ky,s) ~ \^^s, ky, s) + rss'l — kxs'ky, s ) , (12) 

s' 

where eigenstates | — k^g, ky, s) are those from regions and R^. 

The wave function of electrons transmitted to the right reservoir, x > d, is 

s' 

Here r^^' and t^s' are the reflection and transmission matrices. 

The wave functions in the barrier are formed by the eigenfunctions of "oscillating" 
evanescent states (region Bi in Fig. 3): 

lC;^.,.) = E^AA'|Ai^;,A'ir:,fc,), (14) 

AA' 

where K'^ and K'J. are the real and imaginary parts of the wave vector K^, A, A' = ±1. 

Matrices r^g/ , tgs' and b\y are defined by an equation set which follows from the 
boundary conditions [IDl [HI EH] : 

+/?V.^]o_ = Sk' (15) 
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Here the parameter f] = 2Ua/eFz describes the Rashba SOI caused by a lateral electric 
field at the edges of the barrier, is an electric field normal to the 2D layer. The 
SOI constant in the boundary condition disappears because the SOI constants are equal 
all over the sample and the wave functions are continuous at the boundaries. The 
total equation system for the matrices tss', rgs' and b\^, is obtained from the boundary 
conditions for both spin states of incident electrons. Thus, one obtains two systems of 
8 equations each. They are to be added by an equation establishing a relation between 
wave vectors k^^s and K^. This equation follows from the requirement that the energy 
is the same in the reservoirs and the barrier: ({k^g, ky, s) = u + ({K^, ky). 

Now we proceed with the calculation of the charge and the spin currents. Using 
equation f|T3l) . one finds the electron current: 

j{kxs,kyS) = - ^(Xs' - • (16) 

s' 

The spin current is supposed to be defined by the standard expression [161 fT7] : 

J'm = + (17) 

where i = {x, y) designates the current components in the plane, j = (x, y, z) designates 
the spin polarization components, Vi is the electron velocity components. 

The calculation of the spin current in the right reservoir for the state iV'^^j^ys) 
results in the following expressions for the x component: 

where Yf^'^ has following components: 

s' 

X [{Xs + xlKstl.e'^'''-''''^'' + c.c.]/2 , (19) 

s' 

X ts,stls e^^'^--'^-")^' - c.c.]/2 - 2aJ2\tss'\^ 

s' 

- a [tsXAXsX-s + l)e^('=--^-> + c.c] , (20) 

y^x = (kxs + K,)[tsn{Xs + x^)e'^'--'-^^' + c.c.]/2 . (21) 

Here s designates the spin opposite to s. 
The total spin current is: 



s VC^ci 



dk. 



X / , ' j:'y'%CF,ky,s,x). (22) 
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Figure 4. Dependence of the spin polarization of electron current on the distance 
from barrier. The used parameters are: Ep/Eso = 7.99, U/Ego ~ ad ~ "i, (i — Q. 



Straightforward calculations show that the spin current components with polarization 
along X and z directions are absent, Jf^. = J^^ = 0. Only the y component of the spin 
polarization is present in the spin current J^^ ^ 0, just as in the case of SOI absence 
in the reservoirs [10]. The spin current depends on the distance from the barrier. Near 
the barrier, JJ^ oscillates with a period of about 7r/a around a slowly varying value. 
The oscillation amplitude decreases with distance and the spin current asymptotically 
reaches a constant value, as it is shown in Fig. 4. The oscillation is caused by the 
interference of spin-split propagating states whose wave vectors differ by a value of the 
order of a. An electron incident on the barrier with definite spin appears behind the 
barrier in a state which is a superposition of wave functions with different chiralities 
and wave vectors. Their interference results in the spin current oscillations. At large 
distances from the barrier the interference pattern is smeared because the partial spin 
current oscillations lose their coherence due to the dispersion of longitudinal wave vectors 
of incident electrons. The asymptotic behavior of the spin current can be presented 

as 

™ /> X ™ /> X , cosf2aa;+09(Ci7')l 



(23) 



X 

The degree of current spin polarization is quantitatively described by the spin-to 
charge current ratio: 

2 J,%(Cf,oo) 

h JiCp) ■ 

The polarization -P(Cf) calculated as a function of the Fermi energy for two thicknesses 
of the barrier d is presented in Fig. 5. The largest spin polarization is seen to arise when 
the Fermi level lies close to the barrier maximum in an energy interval of the order of 
several Ego- This dependence is similar to that of the case where the SOI is absent in 
reservoirs [10]. This polarization exists at large distances from the barrier. An essential 
point is that due to the oscillation in the vicinity of the barrier, P can by higher or 
lower that the asymptotic value shown in Fig. 5. 

The obtained results depend only weakly on the parameter /?, describing the 
interface SOI. With increasing jS in the range 0-0.1, the general view of the 
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Figure 5. Spin polarization of the current as a function of Fermi energy for two barrier 
thicknesses: ad — 3 (hue 1), ad ^ 8 (line 2). The parameters used are U — 9Eso, (3 — 0- 

P{Cf) dependence remains unchanged, but the degree of spin current polarization 
insignificantly increases. So, P{Cf) increases by 7% as f5 changes from to 0.1. 

4. Spin current transformation in the contact of SOI and normal regions 

In 2D electron gas with SOI, the spin current is known to be a nonconserved quantity 
and therefore its definition is somewhat arbitrary. For this reason an important question 
arises of what quantity is really measurable. In this paper we propose to transfer the spin 
current from the SOI system into a normal 2D electron gas without SOI where the spin 
current is well defined and measurable [Sni EH [22], [23] . To carry out this transformation 
a normal region should be brought into a lateral contact with the SOI system considered 
before. Thus, it is reasonable to extend the discussed system by adding a contact with a 
normal region, which simulates (at least partially) a measuring device. The key problem 
to be solved is to find out how the spin current is transformed while passing through 
this contact. 

The problem is stated as follows. Let a monoenergetic electron fiow is incident from 
the SOI region upon the sharp boundary with a normal region. The spin polarization 
of incident electrons is determined by a non-equilibrium occupancy of spin states at the 
Fermi level, which is characterized by distribution functions of the states with positive 
and negative chiralities, /+(k+) and /_(k_), with k± being the Fermi wave vectors for 
the spin-split subbanbs. One needs to calculate the output spin current in the normal 
region as a function of the spin polarization of the incident current. This problem is 
solved by the same way as in the previous section. Therefore, we describe below the key 
results without going into details. 

Let us consider a simplified case where the distribution functions /+(k+), /_(k_) are 
nonzero only for the states with positive velocity and do not depend on the momentum 
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Figure 6. Dependence of the output spin polarization in the normal region on 

the input polarization in the SOI region in the case of Uc — for different Fermi 
energies (lines 1,2,3). Line 4 shows the result in the case case, when the distribution 
function of incident electrons fills the sector {k^ > 0, < 2a) on Fermi surface. The 
symbols mark the points, in which (/+ — f-)/{f+ + /_) increases from -1.0 to 1.0 by 
steps 0.2, moving from left to right. Inset: the energy diagram of the structure. 



in this sector of the Fermi surface, i.e. /±(k) = f±d{kx). The ratio of spin-subband 
populations determines the degree of the spin polarization of the incident electron flow. 
It is easy to show that 

p(m) ^ TT Cf f+-f~ .24) 

where 

the spin polarization being directed along y axis. 

Note, that far from the contact in the SOI region the spin current does not depend 
on the coordinate x, since electrons occupy the states with well defined spin. Near to the 
contact, but before it, the situation changes essentially because the electrons having been 
reflected from the contact find themselves in a superposition of states with different spin. 
This results in the interference pattern in the spatial distribution of the spin current 
density similar to that shown in Fig. 4. Behind the contact, in the normal region, the 
spin current does not depend on the coordinate since it is a conserved quantity. 

The spin polarization of the output current p(°"*) is defined similarly to equation 
fl23|) as the ratio of the transmitted spin current to the particle current. Of interest is 
the relation between the output polarization, p('"'*)^ and the input one, P^'''\ The input 
polarization is changed by varying the spin-subband population according to equation 
( I24|) . We calculate the output and input polarizations while varying (/+ — /_)/ {f+ + /_) 
to find a dependence of P°"* on P*". This dependence is determined by the Fermi energy 
Ep and the potential step height Uc at the contact between the SOI and normal regions. 
We find that the efficiency of the spin current transformation when transferring through 
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the contact increases with the decrease in Uc- This means that the scattering on the 
contact conributes to the output polarization. The dependencies of p(°"*) on p(*") are 
shown in Fig. 6 for the most favorable case when Uc = 0. They are nearly linear. 
Thus, if Ep ^ Eso, Uc, the spin current passes from SOI region to normal electron gas 
practically without polarization loss. 

Let us address to the problem of the spin polarization of the electron current 
through a barrier studied in the previous section. The largest polarization is reached 
at Ep ~ U. Therefore, if C/ ^ Ego, U^ the spin is transferred into the normal 2D 
gas almost completely even if the distribution functions are uniformly smeared over the 
semi-circle as in the calculation of this section (see lines 1-3 in Fig. 6). In reality, the spin 
transfer efficiency is higher since the distribution function of transmitted electrons /±(k) 
is strongly non-uniform over the azimuthal angle. This occurs because the probability 
for an electron to pass through the potential barrier decreases with increasing \ky\. One 
can say, that electrons are focused by the barrier near to x axis. If the energy is close to 
the barrier top, the characteristic scale of the decrease of /±(k) with ky is of the order 
of 2a. We model this situation by calculating the polarization of the output current 
in the case where the states in sector kx > 0, \ky\ < 2a are only filled. The result is 
presented by fine 4 in Fig. 6. Of course, the scattering processes in the bulk reduce the 
spin polarization because they cause the distribution function to be more isotropic. 

5. Conclusions 

Electron transport through a lateral potential barrier in 2D system with SOI produces 
the considerable spin polarization of the current, with the spin being directed 
perpendicularly to the current. Behind the barrier the outgoing spin current depends 
on the distance in an oscillatory manner, but at sufficiently large distance from it the 
oscillations decay and the spin current reaches a constant value. The most effective 
generation of the spin current occurs when the Fermi energy is close to the top of the 
potential barrier. The maximum degree of polarization at the distance far from the 
barrier exceeds 50%. The spin current generated in the 2D electron gas with SOI can be 
successfully transmitted to a contacting normal 2D electron gas where the spin current 
is unambiguously defined. The spin polarization loss occurring while electrons pass from 
the SOI region to the normal electron gas is neghgible if the contact potential step and 
the characteristic SOI energy are small compared to the barrier height and the Fermi 
energy. 
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